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Abstract: We examine the supersymmetry (SUSY) reach of the CERN LHC operating at 
y/s = 10 and 14 TeV within the framework of the minimal supergravity model (mSUGRA). 
We improve upon previous reach projections by incorporating updated background calcu- 
lations including a variety of 2 — > n Standard Model (SM) processes. We show that SUSY 
discovery is possible even before the detectors are understood well enough to utilize either 
^miss Qr e i ec t rons j n the signal. We evaluate the early SUSY reach of the LHC at y/s = 10 
fNJ ! TeV by examining multi-muon plus > 4 jets, and also lepton-free, acollinear dijet events 

with no missing Et cuts, and show that the greatest reach in terms of ra^ occurs in the 
dijet channel, where it may be possible to probe nig ~ m g ~ 1 TeV with just 1 fb _1 of 
integrated luminosity. The reach in multi-muons is slightly smaller in mi/2, but extends 
to higher values of niQ. We find that an observable multi-muon signal will first appear in 
the opposite-sign dimuon channel, but as the integrated luminosity increases the relatively 
background-free but rate-limited same-sign dimuon, and ultimately the trimuon channel 



yield the highest reach. The optimized reach in these channels extends to m g ~ 600 
(800) GeV for an integrated luminosity of 100 pb _1 (1 fb _1 ). We show characteristic dis- 
tributions in these channels that serve to distinguish the signal from the SM background, 
and also help to corroborate its SUSY origin. We then evaluate the LHC reach in various 
no-lepton and multi-lepton plus jets channels including missing Et cuts for y/s = 10 and 14 
TeV, and plot the reach for integrated luminosities ranging up to 3000 fb _1 at the SLHC. 
For y/s = 10 TeV, the LHC reach extends to nig = 1.9, 2.3, 2.8 and 2.9 TeV for nig ~ m g 
and integrated luminosities of 10, 100, 1000 and 3000 fb _1 , respectively. For y/s = 14 TeV, 
the LHC reach for the same integrated luminosities is to m g = 2.4, 3.1, 3.7 and 4.0 TeV, 
respectively. The reach estimates for ab _1 luminosities may be over-optimistic due to low 
statistics of background with very hard cuts. 
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1. Introduction 



The CERN Large Hadron Collider (LHC) is expected to begin collecting data from pp 
collisions at ~ 10 TeV, with a goal of accumulating 0.1 — 0.2 fb" 1 of usable data in 
the first run. During the very early stages of LHC running (first ~ 0.1 fb -1 ), detector 
commissioning will be in progress, and issues such as detector alignment and calibration 
will be addressed, as the experimental groups use familiar Standard Model (SM) processes 
such as W + jets, Z + jets and ti production to guide the way. Running at 10 TeV is likely 
to continue for a year or more, after which it is expected that the center of mass energy 
will be increased, very likely in several stages, to its design value of 14 TeV. 

While discovery of the Higgs boson (or bosons) or, more generally, the mechanism of 
electroweak symmetry breaking remains a primary goal of LHC experiments, an integrated 
luminosity of ~ 10 fb _1 will be required to claim Higgs discovery (if indeed mHi ggs ~ 
115 — 130 GeV, as indicated by global analyses of electroweak data sets) jy]. An equally 
important objective for LHC is to discover, or exclude, weak scale supersymmetric (SUSY) 
matter. Since production cross sections for strongly interacting sparticles can range up to 
O(10 5 ) fb if rriq ~ m g ~ 400 GeV at yfs = 10 TeV, the hunt for supersymmetric particles 
beyond the reach of LEP2 and Tevatron searches could be very interesting even in the 
earliest stages of LHC running. 

The discovery capability of LHC for SUSY particles is often illustrated with a reach 
plot in the parameter space of some assumed SUSY model || . At each point in SUSY model 
parameter space, many simulated collider events are generated, and compared against SM 
backgrounds with the same experimental signature Q . Judicious cuts are then implemented 
to select out the new physics signals over SM backgrounds, and the signal is deemed 
observable if it satisfies pre-selected criteria for observability. The LHC reach has most 
frequently been analyzed Q within the paradigm minimal supergravity (mSUGRA) model 
||, but other SUSY models have also been studied. 1 

In many models, production of strongly interacting SUSY particles is expected to yield 
the dominant signal channel, at least for rrig ~ rriq ~ 1.7 TeV. Heavy squarks and gluinos 
then decay via complex cascades |IT| which, if i?-parity conservation is assumed, ends in 
the stable (or quasi-stable) lightest SUSY particle (LSP), often assumed to be the lightest 
neutralino Z\. The Z\ escapes experimental detection, so that the generic SUSY signal is 
expected to be the production of multiple high Et jets, multiple high px isolated leptons 
(e or /x, produced via the decays of chargino and neutralino secondaries) and possibly 
also isolated photons, together with missing transverse energy EJp lss . The multiplicity of 
isolated leptons provides a convenient way to classify various SUSY signals [^], and for 
the mSUGRA model, reach contours have been shown for signals in the following channels: 

• jets + E^ 1SS (no isolated leptons), 

1 The LHC reach in anomaly- mediated SUSY breaking (AMSB) models is given in Ref. in the mixed- 
modulus-anomaly mediation case in Ref. jjj], while reach for various model lines in gauge mediated SUSY 
breaking (GMSB) is presented in Ref. |H|, and for gaugino-mediated SUSY breaking in Ref. |]]. The LHC 
reach in mSUGRA with i?-parity violation is presented in Ref. [jl0|. 
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• U + jets + Ef ss , 

• two opposite-sign isolated leptons (OS)+jets + EJp lss , 

• two same-sign isolated leptons (SS)+jets + £™ 1SS , 

• U + jets + £^ iss , 

• 4^ + jets + Ef ss , 

• a real Z -> + jets + £f? iss , 

• a hard, isolated 7 + jets + E™ 1SS - 

These explorations have typically been performed for the design LHC center-of-mass energy 
of y/s = 14 TeV, and integrated luminosities of 10 or 100 fb -1 , anticipated after a year to 
a few years of running at the design luminosity. 

Recently, some attention has been given to the ability of LHC to detect supersym- 
metric matter in the very earliest stages of running when a reliable measurement of E™ 1SS , 
which requires a lead-time for detector alignment, calibration and understanding of the 



performance of essentially all detector components, will not be available [13|. In Ref. [14] 
it was instead suggested that in lieu of E™ ms , high isolated lepton multiplicity could be 
used as a strong cut to reject SM backgrounds at relatively low cost to the expected SUSY 
signal: thus, requiring events with > 4 jets plus OS, or SS, or three isolated leptons could 
allow for probes of rrig ~ 500 — 600 GeV with just 0.1-0.2 fb _1 of integrated luminosity. In 
a follow-up [|l5| paper, the authors refined and restricted their multi-lepton analysis to just 
multi-muons (because reliable electron identification may be difficult in the early stages 
of LHC running). The LHC reach was evaluated for y/s = 10 TeV, and was found to be 
m g ~ 550 GeV for 0.2 fb _1 in the SS dimuon plus > 4 jets channel. 2 Alternatively, a 
search for acollinear dijet events was sugg ested by Randall and Tucker-Smith (RT-S)@ 
as a SUSY search strategy that did not explicitly require E™ ss . By cutting hard on sev- 
eral variables, a signal detectable over SM background was found, especially over portions 
of mSUGRA parameter space where squark pair production is significant and where the 
squarks decay directly into qZ\. 

In this paper, where we re-assess the LHC SUSY reach within the mSUGRA model, 
we have several goals: 

1. We perform much more detailed SM background calculations than many previous 
works, including many 2 — > n subprocesses. In processes such as W, Z and tt 
production, we include exact parton emission matrix elements for the first several 
quark or gluon radiations (see details below). These calculations should model the 
multiple high Et jet production in association with standard processes to much better 
accuracy than the parton shower method. In addition, we include numerous other 
subprocesses such as Zti and tibb production, which have frequently been neglected. 



2 Indeed, muons are already being seen by ATLAS and CMS in cosmic ray events, and further, muons can 
be readily identified at lower pr values than electrons, thus partially compensating for the loss of electron 
channel due to the increased muon signal efficiency. 
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10 ft" 1 


100 fb" 1 


1000 fb" 1 


3000 fb" 1 


y/s = 10 TeV 
y/s~ = 14 TeV 


1.9 TeV 
2.4 TeV 


2.3 TeV 
3.1 TeV 


2.8 TeV 
3.7 TeV 


2.9 TeV (oFIT) 
4.0 TeV (oFIT) 



Table 1: Reach for the gluino mass for integrated luminosity values of 10, 100, 1000 and 3000 fb 1 
at y/s = 10 TeV and 14 TeV, assuming mg - m s . The numbers for 1000 fb" 1 and 3000 fb" 1 should 
really be regarded as upper limits on the SLHC reach. For more details see Sec. [|. 

2. We evaluate the LHC early discovery reach without E™ 1SS cuts in two additional 
multi-muon channels - OS dimuons and trimuons- that have not yet been presented. 
We also show various distributions associated with these quantities that should lead 
to increased confidence that any observed excess arises from a real signal. We also 
evaluate the mSUGRA reach in the RT-S dijet channel, and compare with the multi- 
muon reach. 

3. We present reach plots for the initial energy option of y/s = 10 TeV, and compare 
with similar reach plots for y/s = 14 TeV. We are motivated to do so because just how 
and when the center-of-mass energy of the LHC will be increased to its design value 
is presently unclear: if running full current through the superconducting magnets is 
deemed dangerous, or if it is deemed impractical to re-train the magnets that have 
"lost training" so that these cannot attain the full field, then it may be the case that 
LHC runs below the design energy for the first several years. 

4. We show the LHC SUSY reach for a wide range of integrated luminosities, ranging 
from 0.05 fb" 1 (with cuts pertinent to early reach) up to 3000 fb" 1 with cuts op- 
timized for the extraction of the SUSY signal. These high luminosity values would 
only be accessible at the SLHC, which is intended to upgrade the LHC luminosity to 
C = 10 35 cm s . We stress that several experimental challenges at such high lumi- 
nosities would have to be overcome and our background MC might not be considered 
realistic for such luminosities. 3 Our results for 1000 fb" 1 and 3000 fb" 1 are intended 
to provide an outer limit of the SLHC reach for its first few years of running. With 
this in mind, we show a summary of the LHC/SLHC reach for integrated luminosity 
values of 10, 100, 1000 and 3000 fb" 1 in Table g[ 

The remainder of this paper is organized as follows. In Sec. ^, we present details 
of our improved SM background calculations, along with plots for the total background 
levels at y/s = 10 and 14 TeV from various SM processes. In Sec. ||, we present early 
SUSY discovery reach plots in the various multimuon channels, but with no E™ 1SS cut, 
for several integrated luminosity values and y/s = 10 TeV. We show several distributions 
that would serve to both distinguish the signal from SM backgrounds as well as to make 
a case for its SUSY origin. We also evaluate the LHC SUSY reach in the RT-S dijet 
channel. In Sec. ||, we show updated LHC reach plots for standard mSUGRA signal 
channels including E™ ss cuts and our improved backgrounds, for y/s = 10 and 14 TeV, 

3 For instance, multiple scattering effects will have to be accounted for. 
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Vs = 10 TeV 



Jet and b cuts: 
AR >0.4 
r\ <4.0 

p > 20 GeV 



LO pdf: CTEQ5L 
NLO pdf: CTEQ6L 




Figure 1: Total cross-sections for several SM backgrounds for pp collisions at =10 TeV, with 
three different choices of renormalization and factorization scales (Q) (taken to be equal) shown by 
the solid (Q = V§/6), dashed (Q = V§) and dotted (Q = 2^/1) lines. The NLO results are shown 
as blue crosses. The total cross-section for the SUSY SPSla' and SPT2 mSUGRA cases are also 
shown for comparison purposes. 



and a wide range of integrated luminosities. Because 6-jet tagging, which can potentially 
increase the gluino reach by up to 20% in the mixed bino-higgsino LSP case that occurs in 
the so-called hyperbolic branch/focus point (HB/FP) region of the mSUGRA model |17| , 
will be inefficient at the early stage and is presently uncertain in the ultra-high luminosity 
environment of the SLHC, we do not include it in the present analysis. Our ultimate plots 
include scans over a vast grid of possible cut values, so signal/background is optimized in 
various regions of model parameter space. We end with a summary of our results in Sec. || 



2. Standard model background calculations 

In order to understand how SUSY searches are affected by changes in the beam energy in the 
distinct channels, a careful assessment of the SM backgrounds is necessary. In particular, 
relaxing the E™ 1SS cut may increase the contributions from different background processes 



that are usually neglected in the literature. We used AlpGen[18] and MadGraph[19| to 
compute the following 2 — > n processes: jj, bb, W ± j, Z^j,j^j, bbbb, ti, VV, bbZ, bbW ± , 
bbtt, ttZ, ttW ± , VVV, ttVV, tttt and WW, where j stands for light partons (u, d, s, c 
and g) and V = W^, Z. The leading order (LO) total cross-sections for these processes at 
10 TeV and 14 TeV are shown in Fig. |l] and 0. It is well known that at LO some of these 
cross-sections strongly depend on the choice of the renormalization (^r) and factorization 
(//j?) scales (here we always take fj-F = = Q)- To estimate the systematic error from the 
scale dependence of the cross sections, we calculated these for three different scale choices: 
Q = 2y/s, y/s and \/I/6. As expected, the processes which exhibit a strong dependence on 
the scale are the ones with a cx a™ (n > 2), as seen in Fig. || and ^. In particular, a(bb), 
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mass (CjeVj 


SPSla 


SPT2 


9 


608 


453 


Q 


555 


585 


h 


356 


397 




191 


466 










171 


458 


h 


109 


348 




169 


412 


Wi 


183 


114 




98 


64 



Table 2: Representative sparticle masses for the two mSUGRA case study points labeled SPSla' 
and SPT2 introduced in the text. 

a(tt) and a(bbbb) vary by factors of 1.8, 2.4 and 4.6, respectively. This scale dependence 
is basically the same at 10 and 14 TeV, with a small decrease (~ 10%) for the latter. 
Using MCFM[^] we computed the NLO total cross-sections for tt, Wj, Zj, VV, bbW and 
bbZ, with Q = mt for tt production, and Q 2 = my + Pq^iV) for the the processes. The 
results are shown in Fig. |l| and §. For the dominant backgrounds, namely Wj, Zj and 
ti, the NLO results are well approximated by the LO cross-sections with the scale choice 
Q = vs/6. Hence we choose this scale for all our subsequent background calculations. 
We point out that although our scale choice brings the total LO cross-section closer to the 
NLO result, the same is not necessarily true for the different kinematic distributions used 
in our analysis. However, to be conservative, we do not include a K factor for the signal 
cross-sections. We also show for comparison the total LO sparticle pair production cross 
sections for two mSUGRA points used here as case studies: 

• SPSla': (mo, mi/2, A , tan/?, sign(p)) = (70 GeV, 250 GeV, -300 GeV, 10, +), 

• SPT2: (m , m 1/2 , A , tan/?, sign(p)) = (450 GeV, 170 GeV, GeV, 45, +) 

The SPSla' [^] is a commonly adopted mSUGRA bench-mark point, while the second 
point (labeled SPT2 from here on) has a lighter gluino and slightly heavier squarks than 
SPSla'. 4 Representative sparticle masses for these cases are shown in Table ^. 

Though squark and gluino masses are not hugely disparate for the two points, some 
aspects of the phenomenology are quite different. In the SPSla' case, gluinos decay to 
squarks (with decays to tops and stops occurring about 20% of the time), and qi — > q'W\ 
decays occurring with a canonical branching fraction close to 2/3, and BR(q~R — > qZ\) ~ 
1. For the SPT2 case, squarks mainly decay to gluinos though BR(qi — * q'Wx^) — 
0.3, while gluinos decay via three body modes. The decay patterns of charginos and 

4 The point SPSla' has been selected to yield the correct relic density of neutralino dark matter. The 
point SPT2 has much higher neutralino relic density, but is allowed in scenarios where there exists an 
axion/axino supermultiplet, in which case the DM consists of an axion/axino admixture[^l| rather than 
neutralinos. 
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Figure 2: Total cross-sections for several SM backgrounds for pp collisions at = 14 TeV, with 
three different choices of renormalization and factorization scales (Q) (taken to be equal) shown by 
the solid (Q = V§/6), dashed (Q = V§) and dotted (Q = 2y/§) lines. The NLO results are shown 
as blue crosses. The total cross-section for the SUSY SPSla' and SPT2 mSUGRA cases are also 
shown for comparison purposes. 



neutralinos, however, differ in an important way between the two cases because for the 
SPSla' case, fi, and jj,R are significantly lighter than W\ and Z2, while the sneutrinos 
are just ~ 10 GeV lighter than W\ and Z2. As a result, chargino and neutralino decays to 
stau - remember that the right sleptons, being singlets, have no coupling to winos - are 
significantly enhanced, resulting in a softer spectrum of muons (which frequently result as 
secondaries from r decays) in the SPSla' case. We will see below that this altered cascade 
decay pattern has a significant impact on the early detection of SUSY at the LHC: although 
squark and gluino masses are qualitatively similar, the SPT2 point is accessible at very 
low integrated luminosities, while detection in the SPSla' case requires considerably larger 
integrated luminosity. 

After verifying that for SUSY searches the most relevant backgrounds are tt, Zj, Wj, 
jj, tiZ and bbZ, we improved our results by adding multiple jets to these processes. Using 
AlpGen and the MLM matching algorithm fli~8|| (to avoid double counting) we included in our 
background simulations the following processes: 2, 3, 4 jets, tt + 0, 1, 2 jets, Z + 0, 1, 2, 3 jets, 
W+0, 1, 2, 3, 4 jets, ttZ+0, 1, 2 jets, bbZ+0, 1, 2 jets. In these processes ZW(t*) -> 11, vv (11) 
and — » Iv. Since we apply multijet and hard jet px cuts in our analysis (see below) the 
inclusion of the full matrix element results for the above processes significantly increases 
our background contributions to some of the search channels. 

2.1 Event Simulation 

For the simulation of the background events we use AlpGen and MadGraph to compute 



the hard scattering events and Pythia[23] for the subsequent showering and hadronization. 



The signal events were generated using Isajet 7.78[p3]. A toy detector simulation is then 
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All dimuons (fb) 


(Jb{fj,) (ib) 


bb{fi) (ib) 


tt + jets (Lres) 


60. 7 ± 5.7 


A *7 Q _1_ K A 


0.1 ± 1.1 


Z + jets (Lres) 


OU.O ± (.4 


1 T K J- Q A 
1 ( .0 ± 0.4 


U.U 


XOiai .Dior ^IjrtSSJ 




ft/I ft -r- ft /I 

041.0 zt 0.4t 


K "I _|_ "I "I 


Signal (Lres) 


60.1 ± 0.6 


39.5 ± 0.5 


12.0 ± 0.3 


tt + jets (Dres) 


61.6 ± 6.0 


44.6 ± 5.3 


3.7 ± 0.9 


Z + jets (Dres) 


66.5 ± 6.7 


15.4 ± 3.2 


0.0 


Total BG (Dres) 


128.0 ± 9.0 


60.0 ± 6.2 


3.7 ± 0.9 


Signal (Dres) 


62.7 ± 0.8 


41.0 ± 0.7 


12.2 ± 0.4 



Table 3: Comparison between different calorimeter resolutions for the > 4 jets plus all, OS (with 
a veto for m(^ + /i _ ) < 10 GeV, and 75 GeV < m(iJL + n~) < 105 GeV) and SS dimuon channels for 
the dominant SM background (Z + jets and ti+jets) and the SUSY SPSla' point. The statistical 
(MC) errors are also shown. 



employed with calorimeter cell size Arj x A(f> = 0.05 x 0.05 and —5 < r/ < 5 . The HCAL 
(hadronic calorimetry) energy resolution is taken to be 80%/ i/E + 3% for |?7| < 2.6 and 
FCAL (forward calorimetry) is 100%/\/^+ 5% for \i]\ > 2.6, where the two terms are 
combined in quadrature. The ECAL (electromagnetic calorimetry) energy resolution is 



assumed to be 3%/yE + 0.5%. We use the Isajet[24] jet finding algorithm (cone type) to 
group the hadronic final states into jets. The jets and isolated lepton definitions are as 
follow: 



• Jets are required to have R = \J At/ 2 + Ac/) 2 < 0.4 and Ex(jet) > 25 GeV. 

• Leptons are considered isolated if they have pt(0 > 5 GeV with visible activity 
within a cone of AR < 0.2 of Y,Ef lls < 5 GeV. 

2.2 Hadronic resolution and jet energy scale issues for early discovery 

For our analysis of the early SUSY reach, we considered the possibility that the hadronic 
energy resolution may not be as good as anticipated, which could lead to an underestimate 
of those backgrounds such as Z +jets that fall steeply with Ex{j)- Toward this end, we re- 
evaluated the most important backgrounds to dimuon production, assuming the hadronic 
energy resolution is only half as good as its default value above, i.e. we take, 



• Low resolution (Lres): 160%/V£ + 3% for \q\ < 2.6 and 200%/V^ + 5% for \rj\ > 2.6 
instead of 

• Default resolution (Dres): 80%/VE + 3% for \r]\ < 2.6 and 100%/\/^ + 5% for 
\t]\ > 2.6 

The results are shown in Table [||. 

We see that, with the worse resolution, the Z — > +jets background cross section 

is indeed increased more than the corresponding cross section from the signal, or from 
the top background. However, after the invariant mass cut to veto Z's the difference is 
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no longer striking. We conclude that hadronic calorimetry resolution is unlikely to be an 
issue, even for early detection of a signal. 

One may also be concerned about background uncertainty from the jet energy scale. 
The Z(— > + j cross section is about 100 pb at the LHC, and can be used to establish 

the jet energy scale. The variation of the Z + 4j cross section due to a 5% uncertainty in 
the jet energy scale is ±20% [25|, which yields an estimate of the systematic uncertainty 



for SM background from this source. Since we will require the signal to background ratio 
to exceed 20% for observability (see Sec. ^), uncertainties in the background from the jet 
energy scale also appear to be under control. 

3. Early SUSY discovery: searches at y/s = 10 TeV with no E™ 1SS cuts 

After LHC turn-on in Fall 2009, a period of time will be used for detector studies and 
calibration. During this early phase, scheduled for about eleven months, the LHC will 
operate at yfs = 10 TeV, and accumulate about 100-200 pb -1 of integrated luminosity 



|26, 25]. At this time, the classic SUSY signature of jets + E™ 1SS will almost certainly 
not be viable because of a number of issues related to measurement of missing transverse 
energy E^ 1SS . While weakly interacting neutral particles such as neutrinos or the lightest 
neutralinos that escape detection in the experimental apparatus are the physics origin of 
E™ 1SS , in practice missing transverse energy also arises from a variety of other sources, 
including: 

• energy loss from cracks and un-instrumented regions of the detector, 

• energy loss from dead cells, 

• hot cells in the calorimeter that report an energy deposition even if there isn't one, 

• mis-measurement in the electromagnetic calorimeters, hadronic calorimeters or muon 
detectors and 

• mis- identified cosmic rays in events. 

Thus, in order to have a solid grasp of expected E™ ms from SM background processes, it will 
be necessary to have detailed knowledge of the complete detector performance. Experience 
at the Tevatron suggests that this complicated task may well take some time to complete 
at the LHC because many SM processes will have to be scrutinized in order to properly 
calibrate the detector. For this reason, SUSY searches using the classic jets + E^ 1SB signa- 
ture, or for that matter any signature with i?™ lss as a crucial requirement, may well take 
longer than a year to yield reliable results. 

On the other hand, if sparticles are relatively light - not far beyond the reach of 



Tevatron searches [29] — then their production cross sections at the LHC can be huge, and 
tens of thousands of new physics events may be produced in the first few months of LHC 
operation. For instance, for m g ~ 400 GeV and heavy squarks, the expected gluino pair 
cross sections are in the 10 4 fb range. If trig ~ mg ~ 400 GeV, then SUSY production cross 
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Figure 3: Cross sections for various multiplicities of isolated muons in n-muon +> 4 jet events 
at the LHC, with \/i = 10 TeV. We show the signal levels for the SPT2 sample point by the open 
histogram, along with corresponding levels for various SM backgrounds. In the n(/i) = 2 bin, the 
left, center and right columns show the background components for SS, dimuons and OS (with 
invariant mass cuts), respectively. 



sections are even higher: of order 10 5 fb! Thus, with just 0.1 fb _1 of integrated luminosity, 
we might expect of order 10 3 — 10 4 new physics events to be recorded on tape if the gluino is 
in the 400 GeV range. These large rates provide motivation to re-evaluate search strategies 
that may be reliably carried out at the earliest stages of LHC operation at y/s = 10 TeV. 
To avoid a complicated analysis of the rate at which jets fake electrons which will be rather 
uncertain during early running, we focus on signals involving only muons and jets, and 
where precise determination of the energies does not play a crucial role in the extraction 
of the signal over background. Identification of high pr muons, on the other hand, is one 
of the most straightforward measurements at LHC, and the ATLAS and CMS detectors 
are utilizing cosmic ray muons as a tool for understanding their detectors even before the 
LHC turn-on. 

In what follows, we define the signal to be observable if 

• S > max[5VB, 5, 0.25] 

where S and B are the expected number of signal and background events, respectively. 
The requirement S > 0.25 is imposed to avoid the possibility that a small signal on top 
of a large background could otherwise be regarded as statistically significant, but whose 
viability would require the background level to be known with exquisite precision in order 
to establish a discovery. 

3.1 LHC reach in multi-muon + jets channels without E™ ms requirements 

The center-of-mass energy of 10 TeV is a five-fold increase on the highest collision energies 
currently attained and, as just discussed, represents an opportunity for sparticle searches 
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Figure 4: SUSY reach of the LHC at \/s = 10 TeV via SS-dimuon plus > 4 jets events with 
only the basic cuts detailed in the text, for various integrated luminosities. The fixed mSUGRA 
parameters are Aq = 0, tan/3 = 45 and fi > 0. The solid dots here, and in other subsequent 
figures, denote model points where the signal remains unobservable even for the largest integrated 
luminosity shown in the figure. 



well beyond the reach of the Fermilab Tevatron. Motivated by this, we follow up on earlier 
studies [O, 15 1 and explore the early reach of the LHC in the relatively straightforward 



multi-muon plus multi-jet channels where precise energy measurements are not essential, 
and complications due to jets faking an electron are absent. We impose the following basic 
cuts 5 

• Jet cuts: n(jets) > 4 with E T (ji) > 100 GeV, E T (j) > 50 GeV and \r)(j)\ < 3.0 (jets 
are ordered j\ — j n , from highest to lowest Et) 

• St > 0.2, where St is the transverse sphericity, 

• Muon cuts: p T (fJ-) > 10 GeV, \rj(ji)\ < 2.0, 10 GeV< m(ji + fT) < 75 GeV or 
m(n + fj,~) > 105 GeV (for OS muons only), 

and plot in Fig. || the surviving cross section versus the muon multiplicity for the SUSY 
SPT2 point, along with corresponding contributions from a variety of 2 — > n SM back- 
ground processes at y/s = 10 TeV. At low muon multiplicity, signal is well below the 
background, which is dominated by QCD multi-jet production for = 0, and by ti, W +j 
and QCD production for = 1. For = 2, we see that the signal and background are 
already comparable. We can further divide the dimuon events into the OS and SS class. 
For OS dimuons we apply the invariant mass cuts listed above to avoid the 7*, Z — > 
poles. In this case, the signal (for this sample point) is seen to exceed the background in 



5 Unless stated otherwise, these cuts are imposed on all muon plots in what follows. 
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Figure 5: SUSY reach of the LHC at s/s = 10 TeV via OS-dimuon plus > 4 jets events with 
only the basic cuts detailed in the text, for various integrated luminosities. The fixed mSUGRA 
parameters are Aq = 0, tan (3 = 45 and [i > 0. 



both the OS and SS channels. We also mention that the SS dimuon cross section from 
VF ± ]4 7± + dd/uu production is negligible: specifically, a(pp — > W + W + + X) = 116 fb, 
a(pp — > W~W~ + X) = 46 fb, and this contribution to the + 4j + Aj) cross 

section in Fig. || is just 0.015 fb (0.008 fb). Moving to the 3u. channel, we see that the signal 
drops, but the background, which is dominated by Z, ti and tiZ production, drops even 
further. The 3fi + > 4 jets signal is at the 15 fb level, while the corresponding background 
is around 0.34 fb: despite the fact that the statistical significance as well as the S : B ratio 
are both largest for the trimuon case, an integrated luminosity in excess of 300 pb" 1 is 
necessary to attain the five-event level that we require for observability. 

In Fig. ||, we show the reach of the yfs = 10 TeV LHC for the clean SS dimuon plus 
> 4 jets events for various values of integrated luminosity, using only the basic jet, St and 
muon cuts mentioned above. We scan over the mSUGRA model parameters mo and mi/2, 
with Aq = 0, tan/3 = 45 and \i > 0. We see that with just 0.1 fb -1 , already values of 
nig ~ 450 GeV become accessible. As the integrated luminosity is increased to 0.2 (1) fb -1 , 
the reach increases to 550 (650) GeV. 

In Fig. H], we plot the corresponding reach of the y/s = 10 TeV LHC for OS dimuon 
plus > 4 jets events for various values of integrated luminosity. For low values of m 1 / 2 , 
this signature appears to be even more promising than the SS dimuon channel since this 
signal is observable over portions of the parameter space with m g ~ 450 GeV with just 
50 pb _1 of integrated luminosity! 6 We see, however, that for larger integrated luminosities 
(for which the SS dimuon signal crosses the five-event level), the reach via the SS channel, 

6 We should temper this conclusion with some caution, since it is contingent upon our background 
estimate being correct. In practice, the SM background in the dimuon plus > 4j channel will likely be 
extracted from the data. Assuming that ti events are the dominant source of this background, we estimate 
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Figure 6: SUSY reach of the LHC at y/s = 10 TeV via trimuon plus > 4 jets events with only the 
basic cuts detailed in the text, for various integrated luminosities. The fixed mSUGRA parameters 
are Aq = 0, tan/3 = 45 and /i > 0. 



which has a larger S : B ratio, exceeds that in the OS channel. In the OS dimuon channel, 
the reach in rrig is 500 (600) GeV for an integrated luminosity of 0.2 (1) fb -1 . We should 
remember that the projections for especially 1 fb -1 are conservative, since it is likely that 
by the time this is accumulated, the detectors will be well enough understood for the E™ 1SS 
as well as the electron channels to be useful. 

In Fig. ||, we plot the reach of the -y/s = 10 TeV LHC for the trimuon plus > 4 jets 
events for various values of integrated luminosity. Due to the smaller signal cross-section 
(compared to the dimuon channels), the trimuon signal remains below observability for 
even 200 pb" 1 of integrated luminosity. However, due to the large signal to background 
ratio, even with just 1 fb -1 , this low-rate but relatively background- free channel probes 
gluino masses up to ~ 700 GeV. 

3.2 Characteristics of SUSY multi-muon + jets events 

While a discovery of an excess of SS, OS or 3fi plus jets events would be exciting, it 
would also be useful to check various aspects of these multi-muon events to see if they 
agree with a hypothetical origin from supersymmetry. This is especially crucial for any 
discussion of early SUSY discovery where the signal may initially comprise of just 5-10 
events over a very small SM background. With this in mind, we study various muon (and 



some jet) distributions: any signal in the "counting experiments" of Sec. 3d- will be that 
much more convincing if these events have the expected characteristics discussed below. 
As the accumulated integrated luminosity grows, these same distributions (with electrons 



that an integrated luminosity of 20-30 pb will suffice to extract this background from the measurement 
of the tt cross section with additional jets. 
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Figure 7: A0(/i + ^ ) distribution from OS dimuon + > 4 jets events for SPSla' (dashed) and 
SPT2 (solid) cases, and also for SM backgrounds (shaded). We make no requirement on i?™ lss . 

combined with muons) will, of course, provide more precise information about sparticle 
masses, which will help to zero in on the underlying model. We show these distributions for 
the two SUSY cases SPSla' and SPT2 introduced above as well as for the SM background, 
beginning with OS dimuon channel which has the potential for the earliest discovery of 
SUSY. 

3.2.1 Opposite sign dimuon + jets events 

We begin by illustrating in Fig. ^ the distribution of the transverse plane opening angle 
A(/>(/i + between the muons on OS dimuon + > 4 jets events from the SPSla' (dashed), 
the SPT2 (solid) SUSY cases and from the SM background (shaded). We see that for the 
SPT2 point where a large fraction of the muons originate from high Z<i (— > fJ, fi~ Z\) 
produced in gluino and squark cascade decays; this distribution peaks at small angles. The 
distribution tail comes from dimuons originating from cascade-decay-produced charginos 
and extends out to A0(/i + /i _ ) ~ it. For the SPSla' point, the stau is light, and so a much 
smaller fraction of muons come from direct decays of Z2 , and the corresponding distribution 
is much natter. The difference between the two signal distributions is a reflection of the 
different origins of the muons in the two cases. The SM background is nearly flat, but also 
with a slight peak at low values of A(p. 

The distribution of the invariant dimuon mass m{ji + in OS dimuon + > 4j events 
shown in Fig. |§] is even more distinctive. In this case, the SPT2 signal distribution shows 
the distinctive kinematic mass edge[p7|] at m(/i + /i~) = — = 50.6 GeV. In this case, 
most of the signal dimuons would be expected to cluster just below 50 GeV, strengthening 
the case for the SUSY origin of a signal in the earliest data set. For the SPSla' case 
(which, we emphasize, is somewhat atypical), the mass edge from Z2 — > fi^^ — * ^ H~ Z\ 
decays (which have a branching fraction of just 2.4%) at ~ 82 GeV is considerably less 
distinctive. Moreover, this edge merges right into the Z peak, and so will not be measurable. 
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Figure 8: OS dimuon invariant mass distribution from OS dimuon + > 4 jets events for SPSla' 
(dashed) and SPT2 (solid) cases, and also for SM backgrounds (shaded). In this plot only we do 
not apply the invariant mass cuts for OS dimuons. We make no requirement on E™ lss . 



Furthermore, notice that this distribution is smeared out to low mass values because both 
Z<i and W\ decay dominantly to third generation sleptons, and the muon is frequently 
produced as a secondary from tau decays. Except for the Z-peak, the SM background is 
featureless. 

Finally, in Fig. |9| we show the ^-Br^'ei) distribution for events with OS dimuons 
+ > 4 jets events: the large energy release occurring in sparticle pair production and decay 
provides a harder distribution than that expected from SM background. We see that in 
both cases it should be possible to pick out the signal over the background using only the 
measured transverse energies of the jets, though for the SPSla' case, with the smaller cross 
section, a somewhat larger integrated luminosity will be required. 
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3.2.2 Same sign dimuon + jets events 



As mentioned in Sec. 3.1, the SS dimuon channel requires larger integrated luminosity (due 
to its reduced rate), but provides a much cleaner signal. To get an idea of the overall Pt(m) 
distribution in this channel, we show in Fig. 1C this distribution for the hard (//i) and soft 
{H2) muons from SS dimuon events without any jet cuts for the signal points SPSla' and 
SPT2 and for the background, which mainly comes from ti production. While the highest 
Pt muon from ti production comes from t — > bW followed by W — > jxvu, decay, and is quite 
hard, the lower p? muon must come from b — > cfiis^ decay, and hence is necessarily soft, 
since there is much less energy release in b decays. Thus, while the signal emerges from the 
background only for pt{^\) ~ 100 — 125 GeV, the soft muon from the signal - likely arising 
from some heavy sparticle decay - has a much harder distribution than the corresponding 
background muon. For this reason, we require Pt(P2) > 10 GeV, even though it might 
be possible for LHC detectors to go even lower in muon transverse momentum. We see 
again that the SPSla' case is not the norm in that especially the lower energy signal muon 
frequently arises from tau decay, and so is also soft. 

Next, we show in Fig. [ll] the distribution of the transverse plane opening angle between 
the two muons in SS events with > 4 jets, again for both SPSla' and SPT2 signals and for 
the SM background. This distribution differs sharply from the corresponding distribution 
for OS dimuon events shown in Fig. [7|in that the signal is peaked near A<j) ~ ir. This shape 
is merely a reflection of the fact that in the SS case the two muons typically originate from 
different primary particles in the SUSY 2^2 (gg, gq, or qq) production subprocess, in 
contrast to OS dimuons from neutralino decays. The background distribution is nearly flat 
in A(/>(/i ± /i ± ). 

We have checked that the shape of the distribution of ]P Ex(j) in SS dimuons + > 4j 
events is qualitatively similar to that for OS dimuon events shown in Fig. 0, and so we do 
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Figure 9: Distribution of J^EtU) for OS dimuon + > 4 jets events from signal cases SPSla' 
(dashed) and SPT2 (solid), and from SM backgrounds. We make no requirement on E™ lss - 




p T ((l) (GeV) 

Figure 10: pt distribution for the harder (/zi) and softer (/L/,2) muons in SS dimuon events for the 
SUSY SPSla' and SPT2 mSUGRA cases, along with corresponding SM background distributions. 
We make no requirement on number of jets or E™ lss . 




Figure 11: Distribution of the transverse plane opening angle between the muons, A(/)(/i ± /i ± ), 
from SS dimuon + > 4j events for the SUSY cases SPSla' and SPT2, along with the corresponding 
distribution from SM sources. We make no restriction on E™ lss - 



not show it here. 

3.2.3 Trimuon + jets events 

Because it has an even smaller background, the trimuon channel - at a high enough in- 
tegrated luminosity - could potentially become the best muon channel for SUSY searches 
without E™ 1SS cuts. In Fig. |l^a, we show the trimuon invariant mass distribution from 
3fi+ > 4 jets events for the SUSY points SPSla' and SPT2, along with that from the SM 
background at yfs = 10 TeV. As expected, the signal distribution is relatively featureless, 
since at least two of the muons originate in different parent particles. The trimuon cross 
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section is small - 14.8 fb for SPT2 point, and 3.2 fb for the SPSla' case - so that even 
with 1 fb -1 only the former case leads to an observable signal in Fig. |(| 

The noteworthy thing, however, is that though the signal is very small, it is almost 
background-free, even with just the basic cuts. Thus this channel offers prospects for 
a striking confirmation of new physics (presumably first discovered in the OS dimuon 
channel). In most SUSY models we expect that is comparable to fn^ 2 , so that if 

Z2 can be produced in gluino and squark cascade decays, the chargino can usually also 
be produced via these decay cascades. Thus a subsample of trimuon events is likely to 
include an OS muon pair from Z2 decays. With this in mind, we show the distribution 
of the smaller of the two OS dimuon invariant masses in these trilepton +> 4j events in 
frame b) of this figure. For both cases, we see a mass edge at essentially the same location 
as in Fig. |8|, and further, that below the mass edge, the two distributions are very similar, 
reflecting the common parentage of the dimuons in the two cases. In favorable cases such 
as SPT2, the trimuon signal could thus make a strong case for the SUSY origin of a signal 
first seen with just ~ 50 — 200 pb _1 of the LHC data. For the SPSla' point, the trimuon 
signal offers the possibility of determining the mass edge (obscured by the Z peak in the 
OS dilepton case), though an integrated luminosity of ~ 10 fb -1 may be required (by 
which time it is likely that electron events will also be possible to include in the trilepton 
sample). 

Finally, we remark that the Y2^tU) distribution for 3fi + > 4 jets events is again 
harder than that for the SM background; since its features are again much the same as in 
Fig. 0, we do not show it here. 



3.2.4 Early LHC reach via multimuons plus jets: optimized cuts 



We have seen from the various reach plots presented in Sec. |3,l| that even with very basic 
cuts, and just 0.1 — 1 fb -1 of integrated luminosity, LHC experiments will probe gluino 
and squark masses well beyond the reach of the Tevatron. These cuts have not, however, 



7 Events beyond the mass edge arise, for instance, from g 
Zz decays in the SPT2 case. 



tti decays, as well as from Zs produced in 
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Figure 12: The distribution of (a) trimuon mass, and (6) the smaller of the two OS dimuon 
invariant masses, from trimuon + > 4 jets events at the 10 TeV LHC for SUSY points SPSla' and 
SPT2, along with SM backgrounds. 
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been optimized over different regions of parameter space, and so will not yield the maximal 
reach for a given value of integrated luminosity. In this section, we implement a large grid 
of potential cut values, and then select the grid (set of cuts) value which optimizes the 
reach. We recognize that this may be only of academic interest in that by the time such 
an analysis is actually carried out, LHC detectors may be understood well enough to allow 
the inclusion of electrons as well as E™ 1SS in the analysis. We nevertheless felt that it would 
be worthwhile to explore just how much information can be gleaned from the data, in case 
circumstances make this necessary. 
We begin with a set of pre-cuts: 

• transverse sphericity St > 0.2, 

• for isolated muons: Pt(m) — 10 GeV and \r)(n)\ < 2.0, 

• n(jets) > 2 with E T (j) > 50 GeV and \rj(j)\ < 3.0, 

• for OS dimuons: 10 GeV< m(^+/x~) < 75 GeV or m{yL + pT) > 105 GeV 

Then, for a given point in mSUGRA parameter space with of order 50,000 events 
generated, we find the optimal set of cuts to maximize S/y/B + S, using: 

• n(jets) > 2, 3, 4, 5, 6, 7 



E T {h) > 50, 80 - 340 GeV (in steps of 20 GeV), 400 - 1000 GeV (in steps of 100 
GeV), 

• E T (j 2 ) > 50, 55 - 205 GeV (in steps of 15 GeV), 300, 400, 500 GeV, 

• number of isolated muons n(//) = 2, 3, 4, 5, 6 

We do not implement any .E™ 1SS requirement. We take yfs = 10 TeV, and adopt standard 
parameters Aq = 0, tan/? = 45 and [i > 0. The results of our optimized cuts analysis is 
shown in Fig. ^ for various integrated luminosity choices. We note the following: 

• for low integrated luminosities (0.05 fb _1 and 0.1 fb _1 ) the optimal cuts are n(jets) > 
4, 5 with -Er(ii) ~ 100 GeV and -ErCte) ~ 70 GeV in the total dimuon channel. 

• the same is also true for 0.2 fb _1 and 1 fb , except for ijiq > 1200 GeV and/or 
m i/2 ^ 320 GeV, where the n(jets) > 6, 7 channels become more favorable. Moreover 
the preferred jet Et cuts are always below 200 GeV for 0.05, 0.1 and 0.2 fb -1 and 
below 300 GeV for 1 fb- 1 . 

We caution the reader that the projected reaches will be sensitive to the uncertainty in 
our estimate of the background, especially for high jet multiplicities. It is, however, worth 
bearing in mind that even after optimization, the bulk of the parameter space is probed 



with rij > 4, 5 and = 2. The optimized low ttlq reach from Fig. 13 extends up to 
m i/2 values of 225, 275 and 325 GeV for integrated luminosity values of 0.1, 0.2 and 1 fb -1 
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Figure 13: Reach of LHC for mSUGRA at y/s = 10 TcV for multi-muon +jets events using 
optimized cuts discussed in the text, but without any i?™ lss requirement on the signal. The hxed 
mSUGRA parameters are Aq = 0, tan/3 = 45 and [i > 0. 



respectively. This corresponds to a reach in terms of rrig of about 600, 700 and 800 GeV, 
respectively. 

With accumulation of integrated luminosity, the detectors will rapidly become better 
understood and reliable electron identification will be possible. It will then be possible 
to use different flavor, OS dilepton distributions to statistically subtract chargino and W 
contributions from the same flavor, OS dilepton signal and sharpen up the dilepton mass 
edge first obtained in Fig. ||. Using the same cuts (except that we now include electrons) 
as in Fig. and H we plot the distribution of Acp(e :iz fi^ : ) in Fig. 14, and the "subtracted" 
like-flavor, OS dilepton mass distribution, 



da . , _ , _ , _ _ 
— — (e e + u yfi —eu —e \x 
dm 



in Fig. 03. 



da{e + e ) da(/_i + fi ) 



da(e + fi ) da(e /i + ) 



dm e 



dm 



dm 



fie 



m. 



(3.1) 



3 This has the added advantage that even the "trials factor" is greatly reduced. 
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Figure 14: A0(e ± ^T) for the mSUGRA points SPSla' and SPT2 at ^ = 10, along with SM 



backgrounds. The lepton and jet cuts are as in Fig. [?], and there is no requirement on E, 



miss 
T ■ 



We see that the azimuthal angular distribution of A0(/x ± e =F ) is rather flat, as may be 
expected since the leptons arise from various decay chains including g —* tbW± (SPT2 
case) and g —* tt\ (SPSla' case) which can give e^p^ pairs from the decay of the same 
gluino in addition to e fi^ pairs where the electron and muon each originate in a different 
gluino (or squark) parent (as in the SS dilepton case). As anticipated, the dilepton mass 
edges become significantly sharper upon flavor-subtraction, though the Z peak continues 
to obscure the edge in the SPSla' case. Of course, this subtraction procedure ought to 
work equally well for the trilepton signal shown in Fig. with the understanding that 
the dilepton pairs in (|3.1| ) now refer to the OS dilepton pair in trilepton events with the 
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Figure 15: The subtracted like-flavor, OS dilepton mass distribution for the mSUGRA points 
SPSla' and SPT2 at y^i = 10, along with SM backgrounds. The lepton and jet cuts are the same 
as in Fig. |[ and there is no requirement on E^ lss . 
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Figure 16: Distribution of a = ETij-i) j mijij^) for dijet events with no identified leptons for 
the mSUGRA point m = 350 GeV, m 1/2 = 500 GeV, A = GeV, tan/3 = 45 and /i > 
at ,/s = 10 TeV, along with corresponding distributions from various SM sources. The fixed 
mSUGRA parameters are A — 0, tan/3 = 45 and fx > 0. We require that Ex(ji) + Et{3i) > 700 
GeV, but make no restriction on E^ lss . 



smaller of the two masses. 



3.3 Early LHC reach in dijet channel at y/s = 10 TeV 



Recently, Randall and Tucker-Smith (RT-S)|16| have proposed a search for SUSY in the 
dijet channel, also without using E™ ms . While SUSY dijet + E™ ms searches have been 
proposed for a long time, RT-S emphasized that the search can be made without recourse 
to an E™ ss cut. RT-S propose selecting events with exactly two jets with Et > 50 GeV and 
no isolated leptons. This probes the small mo region of the mSUGRA space because qnqR 
pair production naturally leads to this event topology, since for rriq < m g , q~R mainly decays 
via q~R — > qZ\. RT-S then examine distributions of a) a = Et^H) I m {.3\3'i) ■, b) ^(Pijih) 
(the dijet transverse plane opening angle) and c) the variable £r(ii) + -ErCte)- Signal 
was found to exceed SM background for appropriate intervals of each of these variables, for 
mSUGRA points with light squarks, where the squark-squark and squark-gluino production 
cross-sections are enhanced. 

Following RT-S, we evaluate signal and background for the LHC start-up energy of 
y/s = 10 TeV, and plot the a and Acp distributions for the mSUGRA point itlq = 350 GeV, 
m i/2 = 500 GeV, Aq = GeV, tanj3 = 45 and \x > 0. This mSUGRA point has moderately 
heavy gluinos and squarks (rrig = 1160 GeV, rriq ~ 1000 GeV) and is not accessible via 
early searches in the multimuon channels. After imposing a Er(ii) + Et{j2) > 700 GeV 



cut for both signal and background we see from Figs. 1£ and [L7] that for appropriate cuts 
on a and A</> to remove the enormous QCD background, the signal is above the remaining 
background, which mainly comes from the Z + 2j production, where the Z decays invisibly. 
Dijet events with the required geometry from Wjj and ti production arise only if one (or 
more) of the visible decay products of W or t is missed in the detector. We mention here 
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Figure 17: Distribution of A0(ji, j'2) for dijet events with no identified leptons for the mSUGRA 
point m = 350 GeV, m 1/2 = 500 GeV, A a = GeV, tan/3 = 45 and fi > at y/s = 10 TeV, along 
with corresponding distributions from various SM sources. We require that Ex(ji) + -Er(j2) > 700 
GeV, but make no restriction on E^ lss . 

that for this analysis, we assume that it will be possible to identify the electron in W + 2j 
events if W — » ez^, so that such events can be efficiently vetoed. 9 

To extract the SUSY reach in the dijet channel, we scan over a large range of mSUGRA 
model points and find the optimal set of cuts, using: 

• n(jets) = 2 

• E T (j) > 50 GeV 

• Er(ji) + E T {h) > 100 - 1000 GeV (in steps of 50 GeV) 

• a > 0.05, 0.1 - 0.9 (in steps of 0.1) 

• A(/)(j 1 ,j 2 ) < 0.05,0.3 - 3 (in steps of 0.3) 

• number of isolated leptons (/x or e), n(£) = 

The optimal cuts are selected to maximize S/ y/B + S and satisfy the discovery criteria 
defined at the beginning of Sec. ||. Because the important backgrounds all peak at large 
values of A.<j){j\, J2) , the most effective cut is A</>(ji, j'2) < 2 over most of the parameter 
space. Cutting further on a is then not usually required, so that a > 0.05 is generally 
preferred in order to maximize the signal. As expected, the optimal Et(Ji) + -Et(J2) cut 
increases with m^, but is essentially independent of mo for mo < 450 GeV. 

Our results for the optimized reach in the dijet channel are shown in Fig. [l8|. We see 
that LHC experiments will begin to probe SUSY in this channel even with just 0.05 fb _1 

9 We have also not shown backgrounds from W + j or Z + j production where the W and Z decay to 
taus that decay hadronically. Because of the Er(ji) + Er(j2) cut, these backgrounds will be significant 
only for small values of a, and peak at large values of A<f>, and will be efficiently removed by the same cuts 
that remove the much larger QCD background. 
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Figure 18: Optimized reach of LHC for mSUGRA at y/s = 10 TeV via the RT-S dijet search, 
for various values of integrated luminosity. We assume that it will be possible to veto events with 
electrons or muons, but require no restriction on E™ lss . 



of integrated luminosity, if systematics are under control. With 1 fb , LHC experiments 
should be able to probe to and mo values almost up to 500 GeV, corresponding to 

nig ~ wig just over 1 TeV! In the white region at very small mo, the neutralino is not the 
LSP. As anticipated, the dijet search mainly probes the small mo region where squark pair 
production forms a significant part of the SUSY cross section. 



4. Ultimate reach of the LHC utilizing E™ ss 

As the experiments accumulate data, the detectors will become better understood, and it 
will be possible to utilize both electrons and E™ 1SS (well known to be a powerful discrim- 
inator between SUSY and SM events) in the analysis. Because of technical issues to do 
with re-training of the magnets, at present there is no clear projection for the time-line 
over which the energy of the LHC will be increased from the initial value of 7-10 TeV to 
its design value of 14 TeV. It is clear, however, that the energy increase will be staged 
1 26]. The possibility that the design energy may not be attained for an extended period 



motivated us to study the impact of machine energy as well as integrated luminosity on the 
ultimate SUSY reach of the LHC. Toward this end, we delineate the LHC reach, including 
electrons and E?p 1BB in the analysis of the signal, for y/s = 10 and 14 TeV, for various 
values of integrated luminosity. Projections for intermediate energies may be obtained by 
interpolating between the reach for these two extreme values. 

Before turning to the results, we draw attention to a potentially serious problem that 
arises when we try to make reach projections for integrated luminosities in the ab _1 range. 
In this case, SM backgrounds have to be limited to ab levels to make reliable projections 
for the observability of a signal (near the five event limit) with extremely hard cuts on 
jet and E™ ss . In spite of the very large Monte Carlo background samples that we have 
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generated, we are forced to extrapolate our calculated backgrounds to obtain estimates of 
these backgrounds out to high values of i?™ lss for which our simulation becomes statistics- 



limited. For each process and for each set of jet and lepton cuts listed in Sec. [3.2.4| , we 
extrapolate the background from the lower E?p lss range to higher values of E™ 88 using 
an exponential fit whenever we have enough events in (at least three) low E™ 1SS bins 
to allow a sensible extrapolation. In cases where this is not possible (mostly W and 
Z events with multiple jets and hard jet cuts) we do not extrapolate in this particular 
channel. However to avoid a huge under-estimation of this background when sensible 
extrapolation is not possible, we first check whether the non-extrapolated background 
level (in the channel in question) at low E™ 1SS exceeds the corresponding extrapolated 
background cross-section (often dominantly from ti production which, we have checked, 
can be reliably extrapolated in most cases) by a factor > 5. If it does, we regard this 
set of cuts as unsafe (because the background that we could not extrapolate may indeed 
be too large) and exclude them from our optimization procedure. 10 We recognize that 
this procedure may still underestimate the background at large E™ 1SS in cases where the 
non-extrapolated W and Z backgrounds are below five times the total extrapolated one 
at low E™ ms , but may become the dominant source at large E™ ss , due to a flatter E™ ss 
spectrum (when compared with the extrapolated (ti) one), or simply because we did not 
obtain any W/Z+j events in our simulation. We, therefore, regard the reach obtained after 
applying the procedure just described as the outer limit of the parameter plane that may 
be probed for the corresponding integrated luminosity and label it as oFIT (optimistic fit) 
in the reach plots below. We assume that the QCD background is smaller than the other 
backgrounds for high enough values of E™ ss [28|, so no extrapolation is done in this case. As 
already noted, it is frequently possible to extrapolate the ti background quite reliably, and 
extrapolation of backgrounds from multiple quark and vector boson production processes 
(if needed) is even more straightforward. Finally, we remark that even for these very 
large integrated luminosities, we have moderate control of the extrapolation of the SM 
backgrounds in the SS dilepton and ri£ > 3 lepton channels, where ti production (or for 
higher lepton multiplicities, multiple vector boson processes) is the dominant background 
source and the W and Z backgrounds can be neglected. We regard our projection of the 
reach limited to these channels as conservative, and label it by cFIT in the plots to follow. 

4.1 LHC reach using E™ iss at y/s = 10 TeV 

We begin by reanalyzing the optimized reach of the LHC at yfs = 10 TeV. In addition to 
the optimization over rij, Exiji), Erfo) as i n Sec. 3.2.4 , we now include all leptons (e 



and jj) and also optimize over, 

• E™ 1SS > - 1500 GeV (in steps of 100 GeV). 

• n{&) = 0,1,2,3,4,5,6 

10 We have checked that dropping the set of cuts with the large W/Z background does not affect the reach 
that we obtain in any significant way. The optimization procedure picks out a different configuration where 
the W/Z background is smaller than the top background. 
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Figure 19: The upper frame shows the reach of LHC for mSUGRA at y/s = 10 TeV using optimized 
cuts with N(jet) > 2, electron ID and E™ lss for various values of integrated luminosities. For 
integrated luminosities of 1000 and 3000 fb -1 , the reach is obtained using the oFIT extrapolation 
procedure (discussed in the text) that likely underestimates the SM background, and leads to a 
correspondingly optimistic projection for the reach. The reach in the SS dilepton channel for 1000 
and 3000 fb -1 (labeled cFIT) is shown by the circles in the lower frame, where the corresponding 
oFIT reach is also shown for comparison. We regard this as a conservative estimate of the SLHC 
reach. The fixed mSUGRA parameters are A = 0, tan/3 = 45 and fi > 0. 



where t = [i,e. For the ii£ = 2 OS dilepton signal, we continue to veto events with 
like-flavor, OS dilepton pairs with 75 GeV< m(l+£~) < 105 GeV or m(£ + £-) < 10 GeV. 



Our results are shown in the top frame of Fig. 19. The reach in mj^ for low (high) mo 
extends to = 550 (315) GeV for an integrated luminosity of 1 fb , to be compared 
with 325 (250) GeV in Fig. [If. Thus, ^ iss , together with help from electron ID leads 
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to about a 60% (25%) increase in the gluino mass reach if rriq ~ rrig (rriq ^> rrig). 11 This 
same analysis gives projections for the machine reach with increased integrated luminosity. 
For integrated luminosities of 1, 10 and 100 fb , the gluino mass reach (for rriq ~ rrig) 
is ~ 1.4, 1.9, 2.3 TeV. For integrated luminosities of 1000 and 3000 fb -1 , our projection 
for the outer limit of the LHC reach obtained using optimistic extrapolation of the SM 
backgrounds (oFIT) described above, extends out to 2.8 and 2.9 TeV, respectively. We 
caution that in this region our optimization procedure often selects out the rij = 2, ri£ = 
or 1 topology with very hard jet and E™ ss cuts, where the signal (for the 3000 fb -1 case) is 
just 4 ab. At such low cross-sections the signal is likely termed "observable" only because 
the background from W/Z+j production may be greatly underestimated. Aside from this, 
detector issues that will arise during LHC operation at very high instantaneous luminosity 
have not been included in this analysis, and may significantly reduce the 1000 fb -1 and 
3000 fb -1 reach projections. The oFIT gluino mass reach drops to about 1.3 TeV (1.7 TeV) 
for an integrated luminosity of 100 (1000) fb^ 1 if squarks are very heavy. We have checked 
that, after optimization, the signal almost always arises from gluino and squark production 
except at the largest values of mo where W2 and Z4 production contribute about a third 
of the signal. We repeat that, except possibly at the highest machine luminosities, the 
ultimate reach in m\ji in the HB/FP region will be about 15% higher than shown in this 



12 



figure once 6-jet tagging is utilized to enhance the signal over the SM background [17]. 

To obtain a rough idea of how much the oFIT extrapolation may overestimate the 
reach, we compare the reach - shown by circles and labeled cFIT - that we obtain in the 
SS dilepton channel (for which we have moderate control on the SM backgrounds even 
for ab -1 integrated luminosities) with the oFIT reach in the lower frame of Fig. |l9]. We 
see that for small values of mo, the envelope of even the black circles, roughly speaking, 
follows the 100 fb -1 reach triangles in the upper frame, but for large mo gives a somewhat 
increased reach in miw. Once again, except at the highest values of mo where W2 and 
Z4 production contributes up to a quarter, the signal (after optimization) is dominated by 
squark and gluino production. Modulo detector issues at high luminosity, we regard the 
cFIT reach shown in the figure as a conservative projection of the SLHC reach. 

4.2 LHC reach with Ef iss at y/s = 14 TeV 



We repeat the optimized cuts analysis that led to the reach in Fig. 16, for y/s = 14 



TeV. As for the 10 TeV case just discussed, we first show results using the oFIT to the 



backgrounds for SLHC integrated luminosities in the upper frame of Fig. 20. We see that 



the corresponding gluino mass reach (for rriq ~ rrig), shown in Table [l], extends to rrig = 
(2.4, 3.1, 3.7) TeV for an integrated luminosity of (10, 100, 1000) fb -1 . We see that the 



n We have explicitly checked that the optimization mostly picks out the nt = channel, allowing us to 
conclude that the increased reach is essentially due to the availability of _E™ 1SS rather than of the electron 
signal. 

12 We mention that we have not imposed any cut on the transverse mass between the lepton and E™ 1BB 
in our analysis of single lepton events. Since this cut is very efficient at removing backgrounds from W + j 
and also tt events where just one W decays leptonically, it may be that with an mr (£,E^ iBB ) ~ 100 GeV 
cut, a slightly bigger reach may be obtained in the single lepton channel. In the analysis presented here, 
the 0£ channel almost always yields the largest reach. 
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Figure 20: The upper frame shows the reach of LHC for mSUGRA at y/s = 14 TeV using optimized 
cuts with N(jet) > 2, electron ID and £"™ lss for various values of integrated luminosities. For 
integrated luminosities of 1000 and 3000 fb -1 , the reach is obtained using the oFIT extrapolation 
procedure (discussed in the text) that likely underestimates the SM background, and leads to a 
correspondingly optimistic projection for the reach. The reach in the SS dilepton channel for 1000 
and 3000 fb -1 (and labeled cFIT) is shown by the circles in the lower frame, where the corresponding 
oFIT reach is also shown for comparison. We regard this as a conservative estimate of the SLHC 
reach. The fixed mSUGRA parameters are A = 0, tan/3 = 45 and fi > 0. 



typical gain in the gluino mass reach due to the increased machine energy is a factor 1.3- 
1.4, i.e. this reach, roughly speaking, scales with the center-of-mass energy, and extends 
to over 3 TeV (4 TeV) for an integrated luminosity of 100 (3000) fb -1 . We have checked 
that at the highest m-1/2 an d modest thq values for which we find an observable signal our 
optimization procedure again picks out the rij = 2 channel with very hard jet and E™ 1SS 
cuts with a signal cross section of ~ 2 ab, and for the reasons mentioned in Sec. 4.1, the 
oFIT reach shown is almost certainly an overestimate. 

The most striking new feature in the upper frame is the appearance of a signal in the 
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HB/FP region for an integrated luminosity of 3000 fb _1 . We have checked that the signal 
in this region arises only from W% and Z4 production, in the rij > 2 channel with very hard 
jet and E^ 1BB cuts. Even in the region near mo ~ 4000 GeV and m\/2 = 1200 GeV where 
the SLHC reach contours flatten out, electroweak-ino production accounts for about half 
the signal, whereas for yet lower values of mo, the signal originates essentially in gluino 
and squark events. The signal deep in the HB/FP region is just above the 5a level with a 
cross-section of ~5 ab and, for fixed m^, rapidly falls if mo is reduced because the mass 
gap between the wino parent and the higgsino LSP quickly reduces, softening the E™ 1SS 
spectrum, concomitantly reducing the efficiency for the signal to pass the hard E™ 1SS cut 
required to reduce the background. Indeed, it is very likely that we obtain the reach in the 
HB/FP region only because we get no background from W/Z + j production in one (or 
more) of the many channels examined in the course of the optimization. 

In the lower frame of Fig. we show the corresponding reach in the SS dilepton 
channel. As in the previous figure, the cFIT circles denote points where the signal is 
observable in the SS dilepton channel. Except in the region at m 1( / 2 ~ 1 TeV and tuq > 
3 TeV where about half the signal comes from wino production, the SS dilepton signal 
arises mostly from gluino and squark production. We also see that the HB/FP region 
entirely disappears. Since any signal in this region is most likely to come from W2 and Z4 
production which yields SS dilepton events only if at least one final state lepton fails to be 
identified in the detector, we also checked that there is no observable signal in the ng > 3 
lepton channel, or for that matter in the inclusive SS and rig > 3 event channels. We see 
that as in the = 10 TeV case, our more conservative projection for the reach is again 
close to the projected reach with 100 fb _1 shown in the upper frame of the figure. 

5. Concluding remarks and summary 

LHC experiments will soon probe particle collisions in a qualitatively new energy regime 
and, we hope, uncover new phenomena. In this paper, we access the LHC reach for 
supersymmetric particles, both at the very early stages of LHC operation at the starting 
center-of-mass energy of 10 TeV and limited values of integrated luminosity (0.1-1 fb _1 ) 
when the detectors will not yet be completely understood and calibrated, as well as at later 
stages of LHC operation at yfs = 10 and 14 TeV, again for several values of integrated 
luminosity extending up to the design value, and beyond into the super-LHC stage. To 
obtain these reach projections, we have used improved techniques described in Sec. [2] to 
calculate SM backgrounds from a large number of 2 — > n processes. 

It has, however, often been stated that exploration of new physics (except for reso- 
nances that can be reconstructed as mass bumps) will only be possible after all the detectors 
systems are completely understood, and the "SM is rediscovered at the LHC". This dic- 
tum has been thought to be especially true for the discovery of supersymmetry 13 since the 
experimental determination of E™ ss , which is an essential element of the canonical SUSY 

13 The dictum applies equally well to a variety of proposals that address the mechanism of electroweak 
symmetry breaking and the stabilization of the electroweak scale and include a stable weakly interacting 
massive particle that escape detection in the experimental apparatus fell. 
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signal, truly does require a reliable identification and measurement of all high Ex electrons, 
muons, photons and jets in each event. This motivated us to examine the prospects for 
early sparticle detection without use of E™ ss . 

Since electron fakes from jets could be a serious issue in the early stages of operation, 
and since non-leptonic and single-muon signals will have large SM backgrounds, in Sec. || 
we focused on multi-muon signals from SUSY without any requirement on E™ 1SS . We 
concluded that at the LHC start-up energy y/s = 10 TeV, the OS dimuon channel (which 
has the highest rate of the multi-muon signals) offers the best prospects for the earliest 
detection of SUSY: even with very basic cuts and just 100 pb _1 of integrated luminosity, 
LHC experiments will probe significant portions of mSUGRA parameter space beyond the 
range of Tevatron experiments. The dimuons from the SUSY signal will have masses that 
tend to cluster close to, but below the expected mass edge at — . As more data 
are accumulated, the relatively background-free but rate-limited SS dimuon and trimuon 
channels become the more important, with the latter yielding the highest reach, with 
gluinos as heavy as 700 GeV being accessible with 1 fb _1 of data if nig ~ ms. This reach 
may be extended to 800 GeV using optimized cuts, though by the time such analyses can 
be performed, it is possible that reliable measurement of i? miss will also be available. We 
have also presented a number of muon distributions that serve to characterize the signal, 
and allow us to make at least a circumstantial case for its supersymmetric origin. In this 
connection, see especially the discussion of the low mass dimuon in trimuon events, and 
the distribution of the transverse plane opening angle for OS and SS dimuon events. 

If squarks are light (rriq ~ Trig), LHC experiments will be able to probe squark and 
gluino production in the acollinear dijet channel (mostly from squark pair production with 
squarks decaying directly to the LSP), without the use of £5j? lss |l(| out to masses of several 
hundred GeV with just 200 pb _1 of integrated luminosity, and to ~ 1 TeV with 1 fb _1 . 
We conclude that within the mSUGRA framework, if squarks and gluinos are in the few 
hundred GeV range, there could be a variety of multi-muon and jet signals that will be 
detectable in LHC experiments even during the first run that is expected to accumulate up 
to ~200 pb" 1 of integrated luminosity, despite the fact that the detectors may not be well 
enough understood to allow the use of £™ 1SS or electron signals in these analyses. 

In Sec. ||, we examined the ultimate reach of the LHC for various integrated lumi- 
nosities, assuming that the detectors are fully understood. We show results for y/s = 10 
and 14 TeV in Fig. |l^ and Fig. |2(], respectively. These are succinctly summarized in 
Fig. 21 where, for integrated luminosities of 1000 and 3000 fb^ 1 , we have shown likely 



over-optimistic reach projections obtained using the oFIT extrapolation to estimate SM 
backgrounds with very hard cuts. We observe the following luminosity and energy scal- 
ing rules for the approximate sparticle reach (using optimistic oFIT projections) when 
-mg. 

The reach scales with machine energy, so that the gluino reach increases by about 
40% between y/i = 10 TeV and y/s = U TeV. 

At y/s = 10 (14) TeV, every order of magnitude increase in luminosity gives an 
increase in reach of ~ 400 (~ 600) GeV. We caution though that the oFIT projections 



- 29 - 



N. gts > 2 (with E T cuts, optimized) 





2000 


1 1 1 1 1 1 

3000 fb"' (oFIT) 

"n^OOO^ (oFIT) 


1 i 1 i 1 i 

■' /' s „ 

- ' 

' J' * 

m, =700 GeV . ■•'.<■* 


u 

g 


1500 






✓ — 




10 fl>^ 


\ 




1000 




\ \ 








- ^Sr 

1 ft-' \ 


^-,\3000fb"' (oFIT) 


:: _ - 

lOTeV 




500 


1 


i , i 


14TeV 





I I I I I I I I I I I I I 

1000 2000 3000 4000 5000 6000 

m (GeV) 



Figure 21: The ultimate SUSY reach of LHC within the mSUGRA framework for y/s = 10 TeV 
(solid) and y/s = 14 TeV (dashed) for various values of integrated luminosities. For integrated 
luminosities of 1000 and 3000 fb -1 , we have shown our (over)-optimistic projections obtained using 
the oFIT procedure that tends to underestimate SM backgrounds. A conservative projection for 
the corresponding reach essentially follows the 100 fb -1 contour. The fixed mSUGRA parameters 
are Aq = 0, tan/3 = 45 and fi > 0. Isomass contours for the LSP (double dot-dashed) and for a 114 
GeV light Higgs scalar (dot-dashed) are also shown. The shaded areas are excluded either because 
the neutralino is not the LSP, or electroweak symmetry breaking is not correctly obtained. 



very likely over-estimate the reach for integrated luminosities at the ab level, and 
the growth of the reach with luminosity may well slow down after ~100 fb -1 . 



The situation is more complicated for very large values of mo where signals from 
electroweak W2W2 and W2Z4 production appear to yield an increased reach in the HB/FP 
region at yfs = 14 TeV and ab -1 values of integrated luminosity as seen by the up-turn 
of the black dashed reach contour in Fig. 21. As mentioned earlier, it is very likely that 
the oFIT procedure seriously under-estimates the background from W/Z + j production, 
and that there is really no reach in the HB/FP region. More generally, barring clever 
new strategies (for instance, using gauge boson polarizations) to enhance the new physics 
signal, we believe that the oFIT contours yield over-optimistic projections over most of the 
mo — J7T-1/2 plane. 

Fig. 22 presents a snapshot of the ultimate gluino mass reach of the LHC, with ^/s = 10 
(solid histogram) and 14 TeV (dashed histogram), for several values of integrated luminos- 
ity. The hatched portions of the bar show the LHC reach in terms of rrig. The height of 
the lower hatched portion (bottom-right to top-left hatching) of each bar shows the value 
of rrig up to which discovery of gluinos is guaranteed at the LHC irrespective of the squark 
mass, while the height of the upper hatched part (top-right to bottom-left hatching) bars 
correspond to the maximum value of m g that LHC experiments will be able to probe for 
some value of mo (usually small), where typically rrig ~ m g . For integrated luminosities 
> 1 ab -1 , the reach shown is obtained from the likely over-optimistic oFIT procedure to 
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Figure 22: The ultimate reach of the LHC at </s = 10 TeV (solid) and = 14 TeV (dashed) 
in terms of the gluino and the LSP masses within the mSUGRA framework, for several values of 
the integrated luminosity. Results for the 1000 fb -1 case are obtained using the oFIT procedure to 
extrapolate the background. A conservative analysis would give results close to those for 100 fb -1 . 
The heights of the lower (top-left to bottom-right) hatched bars show the value of nig up to which 
gluinos are guaranteed to be detectable at the LHC regardless of squark masses, while the heights 
of the upper (top-right to bottom-left) hatched bars show the greatest gluino mass that may be 
accessible, usually when to§ ~ m~ g . The dotted bars show the range of the lightest neutralino mass 
over the part of mSUGRA parameter space for which there is an observable SUSY signal at the 
LHC. 



extrapolate the background. Finally, the range of the lightest neutralino mass over the re- 
gion of mSUGRA parameter space where LHC experiments will be able to detect a SUSY 
signal is shown by the red-dotted bars in the figure. Comparing with earlier studies p2| . 
we see that there will be detectable signals in the next round of direct dark matter searches 
- XENON, LUX with ~ 100 kg of noble liquids, or superCDMS (25kg) - over the entire Z\ 
mass range in the figure, and also at IceCube if ~ 550 — 600 GeV, if parameters are 
in the HB/FP region, where the neutralino composition is adjusted to give the measured 
amount of dark matter. Instead, for a bino-like Zi, direct searches with ton-sized detectors 
will be sensitive to Z\ masses up to about 300 GeV. It is exciting that if supersymmetry is 
realized as in the mSUGRA framework (or one of its variants with non-universal parame- 
ters), we expect observable signals not only at the LHC, but also in a completely different 
program of experiments unrelated to accelerator particle physics. 
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